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Abstract. We couple two copies of the supersymmctric mKdV hierarchy by means of the 
algebraic dressing technique. This allows to deduce the whole set of (N, N) supersymmetry 
transformations of the rclativistic sector of the extended mKdV hierarchy and to interpret 
them as fermionic symmetry flows. The construction is based on an extended Ricmann- 
Hilbcrt problem for affinc Kac-Moody supcralgcbras with a half-integer gradation. A gene- 
ralized set of relativistic-like fermionic local current identities is introduced and it is shown 
that the simplest one, corresponding to the lowest isospectral times t±\ provides the super- 
charges generating rigid supersymmetry transformations in 2D supcrspace. The number of 
supercharges is equal to the dimension of the fermionic kernel of a given semisimple element 
E €q which defines both, the physical degrees of freedom and the symmetries of the model. 
The general construction is applied to the N = (1, 1) and N = (2, 2) sinh-Gordon models 
which are worked out in detail. 

Key words: algebraic dressing method; supersymmetry flows; supcrsymmetric affinc Toda 
models 
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1 Introduction 

It is well known that bosonic Toda models are underlined by Lie algebras and that they provide 
some sort of field theoretic realization to them. They are relevant to particle physics because 
they describe integrable perturbations of two-dimensional conformal field theories, allow soliton 
configurations in their spectrum and are useful laboratories to develop new methods relevant to 
the study of non-perturbative aspects of quantum field theory. 

A natural step when having a bosonic field theory is to try to incorporate fermions and to 
construct its supersymmetry extension. In the case of bosonic Toda models this is a not an easy 
task because we want to preserve the integrability, which is one of the main properties of this 
kind of theories. Integrability is a consequence of the existence of an infinite number of bosonic 
Hamiltonians in involution which depend strongly on the Lie algebraic input data defining the 
Toda model itself. Each Hamiltonian generates a bosonic (even) symmetry flow and due to the 
fact that supersymmetry is just a symmetry, it is natural to expect the presence of conserved 
supercharges each one generating its own fermionic (odd) symmetry flow and also to expect 
that the supersymmetric extension is not related to a Lie algebra but to a Lie superalgebra, 
see [12] for an example of how bosonic symmetries are not preserved after supersymmetrization. 
By definition, a supersymmetry is a symmetry where the application of two successive odd 
transformations close into an even one. If there is an infinite number of even flows, then it is 
natural to incorporate the same number of odd flows in order to close the 'flow superalgebra'. 
Hence, the set of fields J- will depend on an infinite number of even and odd variables J- = 
J F (t±i/2,t±i, £±3/2) i±3) • • • )> see [2Q] for a first example of this 'flow approach' applied to the KP 
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hierarchy. Our main motivation to formulate supersymmetric affine Toda models within this 
setting relies on the possibility of using powerful techniques available in the theory of infinite- 
dimensional Lie algebras and integrable systems, in particular, vertex operator representations 
and tau functions. The goal is to set the ground to study the quantization of the affine super 
Toda integrable models within this fashion. 

Several authors have studied the problem of constructing supersymmetric extensions of inte- 
grable hierarchies. On one side, for the Toda lattice most of them use superfields as a natural way 
to supersymmetrize Lax operators while preserving integrability or to obtain a manifestly super- 
symmetric Hamiltonian reduction of super WZNW models, see for example [22 } 1111 \TU[ [5]. The 
common conclusion is that only Lie superalgebras (classical or affine) with a purely fermionic 
simple root system allow supersymmetric integrable extensions, otherwise supersymmetry is bro- 
ken. On the other side, there are several supersymmetric formulations of the Drinfeld-Sokolov 
reduction method for constructing integrable hierarchies in which the algebraic Dressing method 
and the 'flow approach' were gradually developed and worked out in several examples, see for 
example [18^ [9j \19\ H] . The main goal of these works is the construction of an infinite set of 
fermionic non-local symmetry flows but a clear relation between the conserved supercharges and 
its corresponding field component transformations remains obscure. In [3], fermionic fields were 
coupled to the Toda fields in a supersymmetric way in the spirit of generalized Toda models 
coupled to matter fields introduced in [H] and further analyzed in [13] . This coupling was per- 
formed on-shell and only the first half of the supersymmetric sector was analyzed (corresponding 
to the positive part t+i/2, t+i)- An important result of that paper was the introduction of a 're- 
ductive' automorphism r rC( j (constructed explicitly in the sl(2, 1) affine case) devised to remove 
the non- locality of the lowest supersymmetric flow t+i/2; as a consequence, it was shown that it 
is not strictly necessary to start with an affine superalgebra with a purely fermionic simple root 
system in order to get an integrable supersymmetric extension of a bosonic model. See also [23] 
for another (based on Osp(l,4) having one bosonic and one fermionic simple roots) example of 
a Toda model with superconformal symmetry realized non-linearly. The complementary off-shell 
Hamiltonian reduction was developed in [15] by using a two-loop super- WZNW model where 
the (local) action functional leading to the supersymmetric Leznov-Saveliev equations of motion 
was constructed, in principle, for any superalgebra endowed with a half-integer gradation and 
invariant under r re d. It was also shown that several known purely fermionic integrable models 
belong to the family of perturbed WZNW on supercosets where the bosonic part is fully gauged 
away. 

The purpose of this paper is to introduce the second half of the supersymmetric sector 
(corresponding to the negative part t_i/2, i-i) and to study the whole coupled system generated 
by the subset of symmetry flows (t_i, t_i/ 2 , i+i/2^+i)- This analysis was not performed neither 
in [3] nor [15] so this work complement their study. The outcome is that the supersymmetry 
flows described in terms of the algebraic dressing technique turn out to be equivalent to the 
usual notion of supersymmetry described in terms of superspace variables (this is shown by 
considering explicit examples). This allows to locate the supersymmetry of the models inside 
a formalism which is manifestly integrable by construction. 

In Section 12.11 we review the algebraic dressing technique and use it to couple two identical 
copies of the same integrable hierarchy thus defining its extension. In Section [2.21 we introduce 
the relativistic/supersymmetric sector of the extended super-mKdV hierarchy by coupling two 
super-mKdV hierarchies in different gauges. This idea was first used in [21] in the bosonic case. 
In Section 12.31 we construct two odd Lax pairs associated to the chiral sectors of the hierarchy 
and in Section [23] the complete set of extended (N, N) supersymmetry transformations is given. 
The recursion operators are given in Section 12.51 to show that all higher fermionic flows are 
non-local. In Section 12.61 we use the extended Riemann-Hilbert problem to construct a set of 
local fermionic current identities associated to the non-Abelian flows t±i/2i where each pair of 
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isospectral flows t± n is coupled in a relativistic-like manner. It is also shown that the number 
of supersymmetries (i.e. supercharges) is equal to the dimension of the fermionic kernel of the 
operator Ad E. This means that in the superalgebra decomposition g = ker (Ad E) © Im (Ad E) 
induced by the constant semisimple element E, any element of the affine superalgebra g have 
a well defined role, i.e. it defines a symmetry flow or a physical degree of freedom. In Sec- 
tion 12.71 we use the two-Loop super WZNW action to construct the supercharges generating the 
supersymmetry transformations giving a direct relation between the i-ti/2 odd flows and the 
fields transformations. We also show that the Noether procedure reproduces the supercharges 
constructed in Section \2. 61 by using the factorization problem thus confirming their equivalence. 
Finally, in the Sections 13.11 and 13.21 we study in detail the construction in order to have a bet- 
ter feeling of how the fermionic symmetry flows of the models are defined by the kernel part 
(ker (AdE)) and to make contact with the usual notion of superspace. We also give an example 
of a solution to a relativistic-like equation expressed in terms of the higher graded t±3 isospec- 
tral times only, thus generalizing the sine-Gordon equation. In the conclusion we pose the more 
important problems to be treated in the future which are the main motivations of the present 
work. 

2 General analysis 

Here we study the supersymmetric sector of the extended mKdV hierarchy and obtain the main 
results of the paper. The goal of this chapter is to put into one consistent body the new pieces 
with the known previous results. The core of the flow approach we will follow relies on the 
algebraic dressing technique used to unify symmetry flows (isospectral and non-Abelian) of 
integrable hierarchies related to affine Lie algebras. The Riemann-Hilbert factorization defines 
the integrable structure and a related hierarchy of non-linear partial differential equations. 

2.1 The algebraic dressing technique 

+oo 

Consider an affine Lie superalgebra g = g^ half-integer graded by an operator Q 

ieZ/2=-oo 

([Q>0i] = *0i) an d two supergroup elements (dressing matrices) and II taken as the exponen- 
tials of the negative/positive subalgebras of g respectively, i.e. g_ and g+ in the decomposition 
g = g_ +Q+ induced by the projections V±{*) = (*)± along positive and negative grades. They 
are taken to be formal expansions of the form 

e = ex P (v ( ~ 1/2) + + V> ( ~ 3/2) + ■■■), 

U = BM, M = exp ( - V (+1/2) - - V< (+3/2) -•••)> (!) 

where B = expgo € Go an d ip^ £ Qi- The constant semisimple elements E^ 1 ^ of grade ±1 
([Q,E(+V] = ±£( +1 )) define operators AdE^ each one splitting the superalgebra g = + 
into kernel and image subspaces obeying [/C ± ,/C ± ] C K.^, [/C ± ,A4 ± ] C A4^, where = 
ker (Ad E^) and M ± = Im [AdE^ ±l ^). The kernel and image subspaces have bosonic and 
fermionic components /C ± = /C^ © ICp and M.^ = © Mp each one having a well defined 
(half-integer) grade respect the operator Q. 

Recall [6] that the dressing transformation of x € G by g € G is defined by 

9 x = (xgx~ ) , xg~±}. 

The infinitesimal transformation for g = exp A with A = A + + A^ and A± € g± is 

Sax = 9 x - x = ± (xAx~ 1 ) ± x =p xA±. (2) 
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From this we find the pure actions of A = A + G /C+ and A = A^ G /C_ on x = and x = TL, 
respectively 

5 a+ g = -(@a + @- 1 )_@, 5 A _n = + (n^_n~ 1 ) + n. (3) 

To see this, consider A = A + and x = and the upper sign in ([2]). We get 

s A+ e = (eA,.e _1 ) + e - eA + = -(oa+q-^g, 

where we have used the decomposition @A + @~ 1 = (O-A+O -1 ^ + (0 J 4 + @~ 1 )_. For A = A_ 
and x = II the proof is similar. We also have that for A = A_ and x = and for A = A + and 
x = IT the variations vanish, 5 A _Q = and 5 A+ n = respectively. Hence, in the present form, 
the dressing matrices ([1]) only evolve under half of the flows. 

Setting A± = t± n E^ (with [Q, #( ±n )] = ±?ifi( ±n ') and taking the limit t± n -4- 0, we have 
the isospectral evolutions for and II 

d +n e = -(0 J e( +n )0- 1 )_0, d^ n u = +(uE^u- 1 ) + u, (4) 

where 5 A+ @/t+ = (^+0 — 0) /t + — > <9+ n and similar for II. From equations ([4]) we obtain 
the dressing relations 

E^ n) = (OE ( > + ^Q- 1 ) + = GE^O" 1 + <9 +n 00~\ 

E { ~ n) = {jle^tt 1 ) _ = iLE^ir 1 - d^nir 1 

and the Lax operators 

L+n = 9+n — E^~ \ L-n = 9-n + Ejj \ 

The Baker-Akhiczcr wave functions ty± are defined by L± n ^^ = and are given by 
= 0exp I + UnE (+n) J , ^+ = Ilexp ( - ^ t_ n E^ 

\ nGZ+ / \ nGZ+ 

Equations (j4|) describe two identical but decoupled systems of evolution equations as shown 
above, the coupling of the two sectors is achieved by imposing the relation g = \EC ^+ with g 
a constant group element. Alternatively, we have 

exp ( + t+nE (+n) I exp j + £ t_ n E™ ] = Q-\t)Ii{t). (5) 

This is the extended Riemann-Hilbert factorization problem originally used in [2] to extend 
the mKdV hierarchy to the negative flows. From ([5]) we recover (j4]) and two important extra 
equations describing the isospectral evolution of and II with respect opposite flow parameters 

d +n n = +(0£( +n )0- 1 ) + n, d_ n © = -(iiE^ir 1 )©. (6) 

These equations are extended to actions of A + G /C + and A_ G /C_ on II and 0, similar to ([3]) 
we have 

cU + n = +(0A 1 _0- 1 ) + n, ^_© = -(nA_n- 1 )_0. (7) 

The equations ([3]), (jH) and ([6]), (J7]) describe the isospectral evolution and non-Abelian varia- 
tions of the dressing matrices and II and their consistency, as an algebra of flows, is encoded 
in Proposition Q] below. Note that the flows associated to the positive times are dual to the ones 
associated to the negative times, in the sense that K,* + ~ /C_ under the (assumed to exists) non- 
degenerate inner product (*) which provide the orthogonality condition (QiQj) = 5j+j,o of graded 
spaces. This also show how the degrees of freedom are naturally doubled by the extension. 
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Remark 1. If we consider pseudo-differential operators, the equations Q, ([6]) are good starting 
points to extend the KP hierarchy with the negative flows and the expectation value of §5§ to 
extend its corresponding r-function. 

2.2 Relativistic sector of the extended mKdV hierarchy 
From (j3D and (0) we have the following 

Definition 1. The relativistic sector of the extended mKdV hierarchy is defined by the following 
set of evolution equations 

d + e = +(e4 fl) e- 1 ) <0 e, s+n = -(e^e- 1 )^, 

d.e = +(UE { S 1) U- 1 ) <Q e, d-U = -(ue^u- 1 )^, (8) 

for the two isospectral times t±\ = —x^ associated to the grade ±1 constant elements E^ 1 ^ € 0. 
The (*) >0 denote projection onto grades > and the (*) <0 onto grades < —1/2. 

In the definition above we write explicitly the projections (*) ± in terms of grades in order to 
avoid confusion with the different projections used below in (jlOp . The Lax covariant derivative 
(-L = d + A L ^ extracted from ([8]) has a Lax connection A L given by 

L_ = a_ + At A L _ = -B{E { ~ 1) + ^ 1/2) )B~\ 

L + = d + + A L + , A\ = -d+BB- 1 + ^ +1/2) + (9) 

where 

The RHS of ([5]) can be written in an equivalent way because we have 

e _1 (t)n(t) = ®- l (t)BM = (B' 1 e)~ 1 M 

and this motivates the following 

Definition 2. The gauge-equivalent relativistic sector is defined by the following set of evolution 
equations 

d + e' = +(g , ^ +1) g , - 1 ) <0 g', d + w = -(e'£| +1) e'- 1 ) >0 rr, 

d^e' = +(u'E ( s 1) u'~ 1 ) <0 e', d_n' = -(n'Ei _1) n'- 1 ) >0 n', (10) 

where &' = B _1 Q and II' = M. The (*) >0 denote projection onto grades > +1/2 and (*) <0 
onto grades < 0. 

The Lax covariant derivative extracted from (jlOp has a Lax connection 

l'_ = a_ + A't a* = B- X d-B - <pt 1/2) - E ( ~ l) , 

L' + = d + + A%, A+ = B~\E { + 1) + ^ l/2) )B 

and it is related to ([9]) by a gauge transformation L — > L', where A' L = B~ l A L B — dB~ 1 B. 
Clearly, the two definitions are equivalent. 

The constant part of the Lax connection is given by (E is parametrized by x^) 

E^ = E^ 1) dx ± G Q B (E) <8> g (±1) 
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and change under coordinate transformations because of their dx basis. We also have that 

^ 1/2) dx ± = n F (£)®fl (± i/2) 

are fermionic 1-forms. Thus, A L is a superalgebra-valued 1-form. This is to recall that no 
superspace formulation is involved in the construction of our super-Lax operators and that the 
approach relies entirely on pure Lie algebraic properties. 

The equations of motion are defined by the zero curvature of A L , namely [L+, L_] = and 
leads to a system of non- linear differential equations in which the derivatives d± appear mixed 
with the same order, hence the name relativistic. The coupling of one positive and one negative 
higher graded isospectral flow of opposite sign is direct from the construction. This allows the 
construction of relativistic-like integrable equations, see equation (j37j) below for an example. 

In the definitions of the Lax operators above we actually have 

-d + BB- 1 = + Qf, -B- l d-B = AL 0) + QL 0) , (11) 
where (the upper label denoting the Q grading) 

These relations are the solutions to the grade —1 and +1 components of the zero curvature 
relations [L + , L-]_ x = [L' + , L'_\ 1 = for the operators L± and L' ± obtained from ([8]) and (fTUj) . 

The presence of the fermion bilinear results in the non-locality of the odd t±in symmetry 
flows [3] and also in the existence of gauge symmetries of the models as can be deduced from 
the off-shell formulation of the system ([9]) done in [15]. Having Kg ^ translates into the 
existence of flat directions of the Toda potential which takes the models out of the mKdV 
hierarchy. Thus, we impose the vanishing of Q± . Another reason why we impose 

Qi 0) = o, is 

to get a well defined relation between the dressing matrix and the term A^± in the spirit 
of [3], which means that the dynamical fields are described entirely in terms of the image part 
of the algebra M. The kernel part /C is responsible only for the symmetries of the model and 
all this together clarifies the role played by the term Q± . Then, by restricting to superalgebras 
m which = we have local flows and models inside the mKdV hierarchy. 

Remark 2. Flat directions in the Toda potential Vb = {E+ BE_ ^ -B -1 ) allows the existence 
of soliton solutions with Noether charges, e.g. the electrically charged solitons of the complex 
sine-Gordon model which is known to belong to the relativistic sector of the AKNS hierarchy p] 
instead of the mKdV. 

We parametrize the Toda field as B = g exp [ijQ] exp [vC] , provided we have a subalgebra 
solution to the algebraic conditions Q± = 0. The model is then defined on a reduced group 
manifold and (jlip is conveniently parametrized in the image part of the algebra, i.e. 

-d + BB~ l = Af> and -B^d-B = A m . 

The zero curvature (Fl = 0) of Q gives the super symmetric version of the Leznov-Saveliev 
equations [3] 

= Fit 1/2) = -d_vi +1/2) + [B^S 1/S >B-\Et+% 

= F[°l = d^d+BB- 1 ) - [E { + l \BE^ l) B~ l ] - [^ Xf2 \ B^~ X ^ B' 1 ] , 

= F|l 1/2) =B(- d + ^ 1/2) + [E^^'^^B^B- 1 . (12) 
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Written more explicitly in the form 

0- (d + gg- 1 )+d-d + vC = e-*[Ef 1 \gE^g-l}+e-^[^ 
d+ ^J 1/2) = e-^fi^VVi 4 " 172 ^], d-d+nQ = 0, 

the linearized equations of motion with rj = rjo, rjo € M may be written in the Klein-Gordon 
form 

(d+<9_ + m 2 ) o (S) = 0, d+d^v - A = 0, 

for 5 = i// 1 * 11 / 2 ) and log g, where fh? is the mass operator 

m 2 {E) = e- r >°(adEt 1) o adE { + 1] ) o (3) = m 2 IE. 

We have used e -110 [-E 1 ^ -1 ^, E_ ^] = AC. Then, the Higgs-like field rjo sets the mass scale of the 
theory. The massless limit corresponds to t]q —> 00. Note that all fields have the same mass 
which is what we would expect in a super symmetric theory. Taking rj = 770, the free fermion 
equations of motion reads 

d^ 1 ^ = Tm ± (^ Tl/2) ), 

where m ± (*) = e~ v °/ 2 &d E^ 1 ^ o (*). These equations show that fermions of opposite 'chirality' 
are mixed by the mass term and that in the massless limit they decouple. This means that 
positive/negative flows are naturally related to the two chiralities in the field theory. In most 
of the literature, only the positive set of times is usually considered. 

The role of the fields v and r\ associated to the central term C (of the Kac-Moody algebra g) 
and grading operator Q is to restore the conformal symmetry of the models associated to the 
Loop algebra 3 (which are non-conformal) so we are actually dealing with conformal affine Toda 
models. 

2.3 Non-Abelian flows: the odd Lax pairs £±1/2 

Here we deduce the two lowest odd degree fermionic Lax operators giving rise to the ±1/2 
supersymmetry flows, which are the ones we are mainly concerned in the body of the paper. 
The negative part is the novelty here. From ([3]) and ([7]) we have 

Definition 3. The non-Abelian evolution equations of the Dressing matrices are defined by 

s Ki+) e = -(e#)e^) <0 e, s KM u = +(eid+)e- 1 )> n, 

^ ( - ) e = -(n^(-)n^ 1 ) <0 G, s^u = +(mr(-) n - 1 )> n, (13) 

for some positive/negative degree generators and in the kernel of the operators 

Ad £ (±1) . Equivalents, we have 

5^-0 = -(WK^W- 1 )^', s Kl - } W = +(n'id-)n'- 1 ) >0 n'. (14) 

The consistency of all flows, as an algebra, is encoded in the following 
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Proposition 1. The flows (fT3|) and (fTl|) satisfy 

[ S K W> S K (.*)](*) = $\ K W K W](*), [8 K w,8 K l-)](*) = 0, 
where (*) = 0, II, 0', II'. The map 5 : 1C —■ 5/c is a homomorphism. 

Proof. The proof is straightforward after noting that [X, Y]+ = + [-X±,Yy . + 

[X+,Y ± } ± . ■ 

The last relation above means that the symmetries generated by elements in K± commute 
themselves. This can be traced back to be a consequence of the second Lie structure induced on g 
by the action of the dressing group which introduces a classical super r-matrix R=\(T+-TJ) 
defined in terms of the projections V+ and V- of g = g + ± g_ along the positive/negative 
subalgebras g±, see also [9]. The map 8 : K — > 5/c is actually a map (up to a global irrelevant 
sign) to the iJ-bracket (see [7]) [S K ,S K '] = S [K,K'] R , where [K,K'] R = [K,R(K')\ + [R(K),K']. 
Hence, all the symmetries generated by 1C are chiral as a consequence of the second Lie structure. 
In particular, this imply the commutativity of the 2D rigid supersymmetry transformations 
cf. (f20|) below, as expected. 

The ±1/2 flows are generated by the elements ^pD^ 1 / 2 ) £ /C^ 1 ^ of grades ±1/2 in the 
fermionic part of the kernel, where D( ±1//2 ) depend on the infinitesimal constant grassmannian 
parameters. They define the evolution equations (actually variations cf. (J3j) , (|7|)) 

5 +1/2 e = ±(e^ +1 / 2 )e-i) <o 5 5 +1/2 u = -(@d(+^@-% u, 

s_ 1/2 e = -{UD^U' 1 )^, 5_ 1/2 U = +(UD^W 1 )^U, 
giving rise to the dressing expressions 

@{s +1/2 + L>(+ 1 /2)) Q -i = 5+i/2 + D (oj + D (+im = L+i/2) 

U(S_ 1/2 + D^/^U- 1 = 5_ 1/2 + BD^^B' 1 = L_ 1/2 , 

where = [tpi' 1 / 2 ), iK+Va)] € M$ . The derivation of L_ 1/2 follows exactly the same lines 
for the derivation of L+i/ 2 done in [3]. At this point we have four Lax operators £±1/2 and L±\. 
The grade subspace decomposition of the relations [L ±1 / 2 ,L + i] = [-^±1/2? -k-i] = allows to 
take the solution Z?(°) = —5+i/ 2 BB~ 1 . The compatibility of this system of four Lax operators 
provides the 2D supersymmetry transformations among the field components. Indeed, using the 
equations of motion we get their explicit form, see equation ()1T|) and (|18j) below. 
Finally, the odd Lax operators reads 

£+1/2 = K1/2 - 5+i/zBB- 1 + rt +l l 2 \ (15) 
L_ 1/2 = 5_ 1/2 + BD(- 1 MB-\ (16) 

The operator L + i/ 2 was already constructed in [3] and the L_±/ 2 is the novelty here. 

Note that in f)15|) and (|16p are in different gauges. This is the key idea for introducing the 
Toda potential (superpotential) in the supersymmetry transformations which is also responsible 
for coupling the two sectors. 

2.4 Local supersymmetry flows <5±i/2 

The equations (I12p are invariant under a pair of non-Abelian fermionic flows (<5siJSY = ^-1/2 + 
5 + i/ 2 ) as a consequence of the compatibility relations [L ±1/ / 2 ,L + ] = [L ±1 / 2 ,L_] = supple- 
mented by the equations of motion [L+,LJ\ = and the Jacobi identity. They are generated 
by the elements in the fermionic kernel /C^ 1 ^ and are explicitly given by 

5+ 1/2 4t 1/2) = [Et^tB-W+WB], 5+ 1/2 BB-' = 
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<W^ +1/2) = [S +1/2 BB~\ ,/4 +1/2) ] - [d + BB-\ D^)] . (17) 

and 

5_ 1/2 vA" 1/2) = -[B- 1 8_ 1/2 B,^S l/2) ] - [B^B,D^/% 

B-H.yzB = [D^ X ' 2 \^ +1 I% ^ 1/2 VA +V2) = [^.B^-^B" 1 ]. (18) 

The physical degrees of freedom are parametrized by the image part A4 . To guarantee that the 
variations of the fields remain in M we have to check that the kernel components of the above 
transformations vanishes, i.e. 

[5 +1/2 BB~\ Vi +V2) ] G K = 0, [B-^/jjfl, vi" V2) ] G K = 0. (19) 

We will see below in the examples that Q± = imply (|19p as a consequence of the absence 
of the even graded (2n,n £ Z) part of the bosonic kernel Kb in the mKdV hierarchy. These 
conditions turn the lowest odd flows <5±i/ 2 local. 

The Lax operators (fT5|) . (fT6|) generating the odd flows ([17]) . (fTHj) are related to the rigid 2D 
supersymmetry transformations of the type 

N=(N+,N^), 

where N± = dim/C^ 1 ^. As the map L)^ 1 / 2 ) — > <5_|_ 1 ^ 2 obeys 

[<5 ±1/2 , S' ±1/2 ] (*) = (9 T [ ± 1/2iD ,±i /2 ] (*) ~ <9±(*), [<5 +1/2 , <J_i/ 2 ] (*) = 0, (20) 

we see that two fermionic transformations close into derivatives, which is by definition a su- 
persymmetry. This is the case provided \ {F^ ±l / 2 \ F^ 1 / 2 )} ~ E { ± 1] for i^ 1 / 2 ) e /C^ ±1/2) , 
which is significant for the super symmetric structure of the models, see for instance [181. For 
simplicity, we take constant elements E^ 1 ^ which are dual (E^ 1 ^* = E_ , giving rise to iso- 
morphic subspaces /C^ 1//2 ^ ~ /cL 1//2 ^ and to N + = iV_ in consistency with the pairing induced 
by (Ki,Kj) ~ <5i+j. Note that the non-Abelian odd flows close into the isospectral even flows, 
as expected, and that the central and gradation fields do not transform under <5±i/ 2 then, they 
are not truly degrees of freedom of the model. 



2.5 Recursion operators and higher odd flows 

In computing the explicit expressions for odd Lax operators using (fT3|) generating higher degree 
fermionic flows we realize that this is considerably more involved than the ±1/2 cases. Instead 
of that, we use the dressing map /C — > 5fc from the kernel algebra to the flow algebra in order 
to introduce recursion operators. From the relations 

[S K (±l),5p(±l/2)} (*) = <5^(±l) iF (±l/2)](*) = <5p(±3/2)(*), 

we infer the following behavior 

5 F (±n±i/2)(*) = ad^_ (±i) (<5 F(± i/2)) (*) = (TZ ±1 ) n (<5 F( ±i/2)) (*), 
H^i*) = ads K(±1) (*) = fe (± i),*] 

in terms of the recursion operators 72. . The aim is not to reproduce the well known super- 
symmetry transformations but to develop a method to construct systematically all the Higher 
graded odd symmetry flows in terms of its simplest symmetry structure. However, we have to 
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recognize that the use of super pseudo-differential operators and associated scalar Lax operators, 
seems to be more appropriated for computational purposes. 

From this analysis, we have the following two chains of supersymmetry transformations 

X r S K+ c 5 K+ r 5 K+ 

°+l/2 -> °+3/2 0+5/2 ->■ °+7/2 ~> ' ' ' > 

■■■ *~ 0-7/2 $-5/2 £-3/2 0_ 1/2 , (21) 

where the ones corresponding to <5±i/2 are considered as starting points. The variations 5j^± are 
given by (|13p or (|14p . For example, for a degree +1 element K^ +1 ^ we have from f)13|) that 

<W)^T 1/2) = -[^- 1 \Br\Q^Q-% 1/ JB\J. (22) 

The dressing matrix G factorizes as Q = C/S 1 , where £7 £ exp.M is local and S 1 € exp/C is 
non-local in the fields [3], splitting the Dressing of the vacuum Lax operators (L± = GL^G" 1 ) 
as a two step process. A U and an S rotation given respectively by 

U^L+U = d + + E ( + 1] +K { ~\ XJ- X L-XJ = 8- + K { ~\ (23) 
S- 1 (8 + + E { + 1] + K^) S = d + + , S- 1 (cL + K^) S = cL + £ M) , (24) 

where K± ' G IC involve expansions on the negative grades only. The components tp^\ i = 
— 1/2, —3/2, ... of U are extracted by projecting (j23[) along AA and the components sW, i = 
— 1/2, — 3/2, . . . by projecting (|24j) along IC. This allows to compute (|22p . The higher graded 
supersymmetry transformations are inevitably non-local because of the presence of the kernel 
part S appearing in the definition of the transformations S K (±i) used to construct them. Thus, 
the best we can do is to restrict ourselves to a reduced manifold (defined by Q± = 0) in which 
6±i/2 are local. From (fT3]) we have 



[S K (+i),6_ 1/2 ](*) = <y[ A -(+i) ) _ D (-i/ 2 )] (*) = 



and we cannot connect o~_i/ 2 and 5+1/2 through a flow, reflecting the chiral independence 

of the <5+i/2 transformations as a consequence of the i?-bracket. This is why in (|2ip the sectors 
are treated separately. Although the higher graded odd flows are non-local, their square always 
give a local even flow. A similar conclusion for this behavior was found in [8] by using superspace 
formalism. 

2.6 Generalized relativistic-like current identities 

In this section we derive an infinite set of identities associated to the flows generated by /C^ 1 ^ . 
The word relativistic is used in the sense that each t± n is coupled to its opposite counterpart t Tn . 

Proposition 2. The infinite set of fermionic local currents defined by 

4+V2) = (£»(+V2) es (+n) e -i) ) j(+V2) = ( jD (+i/2) n ^(-)n- 1 >, 

J™ = (£(-l/2) e /£(+n) Q /-l^ jM/2) = ^ D (-l/2) U , E (-n) uf -l^ £5) 

satisfy the following identities 

d +n &* - d. n J^ = 0. (26) 
The e Kl ±1/2) are the generators of the fermionic kernel. 
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Proof. The proof is extremely simple and is based only on the relations Q and ©. Start with 
d+nJ^ = < J D (+1/2) [(G J E(+™)e" 1 ) >0 , (ILE^II- 1 ) 



d-nJ { +^ /2) = ( J D (+1/2) [(ei? (+m) e- 1 ) >0 , (rLE^n- 1 ). 



'<0J /' 

l<J> 

to get 

- < j d(+ 1 /2) [(QEi+^e- 1 ) > , (m^rr 1 ) <0 ] >. 

This sum vanishes for m = n. For j( _1//2 ) the proof is analogous. ■ 

These identities mixes the two sectors corresponding to positive and negative isospectral 
times in a relativistic manner. They can be written in a covariant form r] ij ' J^- J^ 1 / 2 ' = if we 
define a constant 'metric' rj = rjijdtidtj for each pair of positive/negative times. However, the 
interpretation of this higher graded 'light-cone coordinates' deserves further study. 

Consider now the lowest isospectral flows t±i = —x^. The current components (|25p are given 

by 

ji+l/2) = _( £ ,(+l/2)[^(-l/2) jS+SB -l]) | jC+1/2) =+ ( D (+M B fJ; ( L im B- 1 ), 
J™ = +<D(-V2) B -1^(+V2) B>) j(-1/2) = + ^(-l/2)^(+l/2) >S -l a _ S ]). 

Then, there are N = dim/C^ tl//2 ' ) associated relativistic conservation laws (for each sector) given 
by d + J^ ±:L / 2 ) — cL J^ 1 ' 2 ' = 0. More explicitly, we have 

d-{l^V 2 \d + BB^]\ K ) +d + {B^- 1/2) B-\) =0, 

-a_(S-Vi +1/2) B| K ) +a + (^ (+1/2) ,i?- 1 5_ J B]| x; ) =0. (27) 

This time, the identities provide supercharge conservation laws due to the fact that the flows t±\ 
are identified with the light-cone coordinates x^ = ^ (x° ± x 1 ). It is not clear if the identities 
associated to the higher flows t± n , n > +1 provide new conserved quantities because one is not 
supposed to impose boundary conditions or integrate along these directions. For higher times 
they are taken as simple identities consequence of the flow relations above. 

Now that we have N = dim/C^ 1 ^ supercurrents associated to /C^ 1 ^, let's compute their 
corresponding supercharges by the Noether procedure in order to check that they really generate 
the supersymmetry transformations ()1T|) and ()18|) . 

2.7 Supercharges for the SUSY flows <5±i/2 

The action for the affine supersymmetric Toda models was deduced in [15] and it is given by 



STodat^] = S wzmv [B] - A f <^ +1/2) 5_V ( - 1/2) +^ 1/2) 9 + ^ (+1/2) > 

+ A f (E^BE^B- 1 + ^ 1I2) B^ 1I2) B- 1 ). (28) 



This corresponds to the situation when we restrict to the sub-superalgebras solving the condition 
= 0. In this case the potential ends at the second term providing a Yukawa-type term turning 
the model integrable and supersymmetric. The light-cone notation used for the flat Minkowski 
space £ is x = \ (x° ± x 1 ), d± = do ± d±, rj^ = rj |_ = 2, j/" 1 = r\ h = i e_| = — e |_ = 2, 
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e + = — e + =2 corresponding to the metric 7700 = 1, f]u = — 1 an d antisymmetric symbol 
eio = — eoi = +1. A coupling constant is introduced by setting E^ 1 ^ — > fJ*E± and ^( ±1 / 2 ) — y 

/i-V2^(±l/2). 

An arbitrary variation of the action (|28|) is given by 

^S Toda = / (SBB^Ffl)- [ (S^+^B-^y^B)- [ <^(-i/2) F CH/ 2 )> 



and the equations of motions are exactly the super Leznov-Saveliev equations, cf. (|12|) above. 

Taking 5 — >■ <5susy = £-1/2 + ^+1/2; using (fTT|) . ([18]) and considering Z)( ±1 / 2 ) as functions of 
the coordinates a; , we have the supersymmetric variation of the action 



2vr 

T 



W^Toda = / s (d-D^iB-^^B) - d + Dl-V*> [^ 2 \B-^B]) 

-J (d-D^w [^(- i /2) )5+jBS -i] +a +J D( +i / 2 )(i?vA" 1/2) i?- 1 ); 



This allows to obtain two conservation laws 

= / (D^(dj + ±i/2) +d + jy /2) )), 

which are exactly the ones derived by using the extended Riemann-Hilbert approach (|26p for 
the lowest flows (|27p . Then, there are dim/C^ supercurrents and supercharges given by 
flow 5+1/2: 

A- 1/2) = y^ 2 \d +BB yy ym = ^( 7 V2) jB -i kj 

q+ = J dx i {[4- l/2 \d + BB- i ] + B^y /2) Byy (29) 

flow 8-\/2'- 

A +1/2) = -bt^wb\„ y +l/2) = [4 +1 ' 2 \b-^b] y 

q- = J dx i {y +i i 2 \B- i d^B]-B-y { y /2) B)y (30) 



The variation above is the same when (|19p are zero or not, this is because all the fields 
are defined in A4 and the kernel part does not affect the variation at all. These two ways of 
extracting the supercharges show a deep relation between the algebraic dressing formalism and 
the Hamiltonian reduction giving (|28D . 

Now specialize the construction done above to the simplest toy examples. The supercharges 
are computed from the general formulas (|29p and (|30p . We want to emphasize that the sub- 



superalgebras solving the condition = have no bosonic kernel K-b of degree zero in 

consistency with the absence of positive even isospectral flows t+2n. m the mKdV hierarchy. 



3 Examples 

These examples show how the superspace notion of supersymmetry can be embedded consistently 
into the infinite-dimensional flow approach. The usual SUSY transformations corresponds to 
the flow algebra spanned by the times (t_i, t_i/2, £4-1/2 > £4-1) • We can have several pairs of odd 
times t±i/2 depending on the dimension of JC^ 1 ^ 2 ^ as shown above. 



Supersymmetry Flows of Affine Toda Models 



13 



3.1 The N = (1, 1) sinh-Gordon model reloaded 

'[1] 



Take the si (2, 1)L superalgebra (see Appendix [A] for details). The Lagrangian is 



L = -— [d+td-t + + - V] , 

V = cosh [2(f)] + 4^+ V>- cosh [(f)] (31) 

and the equations of motion are 

3 + d-(/)= -2n 2 sinh[2<£] - 2/i^+^- sinh[0], 

= 2/^/>_ cosh[</>], <9 + V'- = —2^ + cosh [(/>]. (32) 

With _D( +1 / 2 ) = e-Fg 1 j = e+F^ 1 ^ 2 ' and tp± — > hif)±, the supersymmetry flows 

are 

8±i/2(t> = ±e T ^ ± , <*±i/2^± = Te T d ± 4>, S ±1/2 ip T = 2fie T sinh [<?!>], 
where we have used the parametrizations 

B = exp^], V (+1/2) = ^G{ +1/2 \ V ( " 1/2) = i>+G { i 1/2 \ 
^~ l/2) = 2^G 2 " 1/2 \ V + +1/2) = -2V+Gi +1/2) . 



We can check (|20p by applying the variations twice giving 

[*±i/ 2 ,*±i/ 2 ] =2e T e' T a ± , [^+1/2^-1/2] =0- 

Then, we have two real supercharges N = (1, 1) because of dim/4 ±1/2) = 1. They arc 
given by 

*±i/2 : Q iTl/2) = Q±F[f /2 \ Q± = J dx 1 (i;±d ± (j) T ^ T ti(<f))) , 

where h{4>) = 2// cosh [0] and h'(4>) is its functional derivative respect (f>. 

Now rotate the fermions by the phase exp (in /4) in order to write (|3ip in a more familiar 
form 



z % 
-d + (f)d-<p + -ip-d+ip- + -V'+<9-?/'+ - V 



L = -* 

7T 

V = ^{ti(q>)) 2 +ih"( ( f ) ^ + iP_+ f i 2 , 



which is known to be invariant under the iV = (1, 1) superspace transformations for a real 
bosonic superfield. The area term comes from squaring h'((f)) = 2/zsinh[</>]. 

Note 1. The Poisson brackets are defined by 

where e = 1, for bosonic- fermionic quantities and vrj = gnjfjs ■ The Dirac bracket is defined by 

{A, B} DB = {A, B} PB - {A, &} PB (C- 1 )^ B} pB , 
where Qj </>i} PB an d 4>i are the second class constraints. 



14 



D.M. Schmidtt 



With the Dirac brackets {(j),dt4>} = 1, {V'iiV'i} = ~~ * an d Q±f = {Q±, /} we have, after 
replacing Q+ — > iQ±, the action of the supercharges on the field components 

Q±4> = -iil>±, Qhl>± = +9±<A, Q±i>T = Tti(4>). 

Finally, the total flow can be written as <5susy = — ie~Q\ + ie + Ql_ 
5svsy4> = +ie-ip + - ie+</>_, ^susvV'i = T£^d±<t> + e±h'(<p), 

which are the ordinary N = (1,1) supersymmetry transformations obtained by using a super- 
space approach. 

Now, construct the £+3/2 transformations starting from £4-1/2 by applying 5 K (+i) as shown 
in (|2ip . From [3] we have 

$k(+i)<I> = -^ip + Q(x + ), S K (+i)if} + = +-a (<9+V>+ - d + (j)Q{x + )) , 

S Ki+1) ^ = -la (ip+htf) - h'(4>)Q(x + )) , 

(V>+<9+0). From the relation [K { 2 +1 \ F 2 (+1/2) ] = F x (+3/2) we have that 
£+3/2 = [^r-(+i)j ^+1/2] an d the transformations are given by 

5+3/2^ = e- " +( f>Q(x + ) + l^ + H{x + )^j , 

5+3/2^+ = e- (ld 2 +(f> + d + ^ + Q{x + ) - ld + <pH(x + )^j , 

<5+ 3/2 0- = e_ (-±0+^(0) - i/>+Q(x + )h(4>) + ^'(0)^(x + )J , (33) 

where lf(x + ) = ^ (^(d + <j)) 2 + ^+(9+V>+) • We also find the variations 
Ki/2Q{* + ) = ~e-H(x + ), 

S + 3/2Q(x + ) = 6- (I (d + <P) 2 - l^+d+i; + + d + ^ + Q{x + ) - ^ 2 (* + )) , 

5 +3/2 H(x + ) = e_ (d + 4>d + ^ + - ld 2 + W+ + ld + ^ + H{x + ) - Q(x + )d + H(x + )^j . 

Applying £+3/2 twice we get a local flow description of the hierarchy for £+1 and £+3 in terms 
of the sinh-Gordon variables used to described it in terms of t+i and cf. equation (|32p 



49 +3 <A = ^ - 2 (8 + <pf - 3d+H>+d+il>+, 4cW+ = d%^+ - 3<9+</><9+ (d+H>+) , 
49+3 tp- = (2 (O+0) 2 - d 2 + ^+) , (34) 

where [£ +3/2 , <5 +3/2 ] (*) = -2eie 2 d+ 3 (*) in agreement with {F$ +3/2 \ f[+ Z/2) } = -2E { + 3) . In- 
troducing u = d+4> (the space variable x is described by i+i) we recover the super mKdV 
equations 

4d +3 u = u'" - 6u 2 u' - 30+ (uip' + )' , 49+3-0+ = ^+ - 3u (uip + )' , 

4<9+ 3 V>- = (2-u 2 0+ - *l>+) Kd^u). 

The V- equation is non-local and is remnant from the negative part of the hierarchy. The 
supersymmetry has to be reduced to the usual N = (1, 0) in order to have a local description in 
terms of the mKdV variables u and if)+. Note that Q{x + ) and H(x + ) resembles a piece of the 
supercharge and a component of the stress tensor. 
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Remark 3. The higher grade fermionic transformation are some sort of 'square root' of the 
isospectral flows, this can be seen also in jS] where Poisson brackets were used (compare 
with (jUD ). 

The positive and negative parts of the extended mKdV hierarchy carries exactly the same 
information when considered separately and to obtain relativistic equations we combine their 
Lax pairs in different gauges as explained above by general arguments. The potential couples 
the two sectors and enters through the functional h((j>). We check this explicitly by considering 
(i+i, £+3), (t-i, i-3) and compute L± 3 in order to construct a 1-soliton solution which solves any 
equation of the hierarchy mixing the four times (i+i, £+3, t-i, £-3)- We keep the name soliton 
solution for simplicity but strictly speaking we need to have a multi-vacuum theory in order to 
define the asymptotic boundary values of the field solution. 

The algebraic dressing technique (recall equation Q and (J6j) ) suggest the following forms for 
the Lax operators 

L +3 = d +3 + + + Z)( +1 ) + + D (+2) + £,(+5/2) + E {+S) i 

= d_ 3 - B{D^W + + D i-m + D (-2) + £,(-5/2) + E i-^ B -\ (35) 

From [L + ,L + 3] = we get for the positive part of the hierarchy, following [3], the solution 

£.(+5/2) = (-2^ + )GS +5/2) , = (-d + <t>)M[ +2 \ 

= (-2^+8+^+) K[ +1) + Q {d + 4>f + 2^ + a + V+) + Q«9 2 A M 2 (+1) , 



.(+1/2) 



£(0) 



The equation of motion are given by the degree zero component and are, as expected, given the 
first two equations of (|34p . after taking — > Now, performing a gauge transformation 

with I? -1 in order to eliminate the B conjugation on L_3 and L_ we change to 

L_ = a_ + S^cLS - (£ M) + ^ 1/2) ) , 

£_ 3 = 5_ 3 + £-i 5 _ 3 i? - {D^W + Z)(- 1 ) + Z)(- 3 / 2 ) + D(- 2 ) + J D(~ 5 / 2 ) + £(~ 3 )) 
and from L_ 3 ] = we get for the negative part of the hierarchy, the solutions 
= (2^)G 2 ~ 5/2) , = {-d^)M[~ 2 \ 

= w-d-fi f 2 ( " 3/2) + (a_^_) gS" 3/2) , 

DC- 1 ) = (-2</>_d_^-) + Q (9-0) 2 + 2^_a_^_J + (—d-A Mt'\ 

D(- 1 /2) = _1 (q2_ h _ _ Q_ ( j ) Q_ 1 p_^ F ("V2) + Qd 2 V- " (9-<A) 2 G 2 " 1/2) - 

The equations of motion are now 

4(9_ 3 </> = &L<t> - 2 (d_<^>) 3 - Sd-fap-d-ip-, 4d_ 3 V>- = - 35_0S_ {d-H)-) 

after taking ?/>_ — > With u = c?_0 (the space variable x is described by we get 

Ad- 3 v = v'" - 6v 2 v - 3V>- (W*-)' , 4a_ 3 ^- = ^_ - 3t> («^_)' . 
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From the solution of the equation (03]), i.e. = —d+ 3 BB~ l we confirm that 49+30 = 

d\cf) — 2 (d + (j)) 3 — 3d + (pip + d + ip + , after taking ip + — > \^+- This is the N = (0, 1) mKdV equation 
with opposite chirality. 

The simplest example of a relativistic equation is provided by the sinh-Gordon model which 
corresponds to the lowest t±± times. Now compute [L +3 , £-3] = with L_3 in the form (f35|) . i.e. 
in its original gauge because the two copies are to be taken in different gauges. The generalized 
relativistic system of equations is given in the Appendix|Bl equation (|44p . and is a generalization 
to the £±3 times. The two sectors of the extended hierarchy are identical, thus we complete the 
system (|34p with the following set of equations 

4d_ 3 </> = d 3 _(/)-2 (d_0) 3 - 3d_c^-d_V>-, 4<9_ 3 V'- = dLrp- - 3<9_c/>d_ (<9_<^_) , 
45-3^4- = ( 2 {d-4>? V>- ~ d^%b-)h((f>) (36) 

and due to this symmetric behavior we can easily read off the £-3/2 transformations directly 
from (|33p simply by replacing all + sub-indexes by — sub-indexes. 

Now consider for simplicity the bosonic limit of the generalized equations (j44]). It is given by 

d +3 d_ 3 = -2sinh[20] + -( (d + <j)) 2 {d^f + d%(j)d 2 _cj)) sinh[20] 

- ^(dl^id^) 2 + (d + <f>) 2 d 2 _cj>) cosh[20], 

d±z{d T cp) = (d ± 0) 2 sinh[20] - d 2 ^ cosh [20], 

d±3 (d T (/)) 2 = 2<9 T 0(d ± 0) 2 sinh[20] - 2d T 4>d%4> cosh [20], 

d±3 (d^cfi) = 4d ± + 2d T 4> (d±(j)) 2 cosh[20] - 2<9 T 0d|0sinh[20] (37) 

and describe the behavior of the '(^-descendants' <9±0, (<9±0) 2 and d±cft in terms of opposite £+3 
times. The equation (|37|) can be obtained alternatively by using the basic non-linear relations 
d±3(f> = \d±4> ~ \ (9±<i>) 3 , and the lowest relativistic equation <9 + <9_0 = — 2sinh[20]. The whole 
set of equations associated to the times (i_3, t+i, ^+3) is completed by the equations ex- 
tracted from the relations [L + ,L^ 3 ] = and [L^,L +3 ] = 0. They are given by the sinh-Gordon 
equation itself, the equations on the third line of (|37p and 

43 T = d± {d T <j)f sinh[20] - <9±d 2 0cosh[20], 

= d±d 2 , 0sinh[20] - d± (d T 0) 2 cosh[20]. (38) 

Now take K ±{ ^ = E^ in flT3]) and conjugate them with the grading operator Q as follows 

£'(±") = exp(aQ)£( ±n ) exp(-aQ) = exp(±n)£( ±n ) = A ±n E^ ±n \ 

where A = exp(a). The equations are invariant under these rescalings and the Lorentz transfor- 
mations, i.e. x = A^x^ can be generalized to the whole set of flows by taking t± n = A ±n t± n . 

To find the 1-soliton solution of the equation (|37p we use the dressing method which is 
another application of the formalism given above in Section 12.11 see also |16] for computational 
details. The four vacuum Lax operators involved are 

L± = d± ± E^ , L ±3 = d± 3 ± E^ 

and the zero curvature conditions imply they are pure gauge AY = T^diTo. Hence we have as 
usual 

Al = - x-C, A v _ = -E^\ Alt = £(+ 3 ) - 3t_ 3 C, AY_ 3 = -E^ 3 \ 
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T = exp (i +3 £ (+3) + exp ( - x~ E^ - t^E^), 

where we have used \E^ m \ E^] = \{m — n)5 m +n,oC for (m, n) odd integers. The dressing of 
a vacuum Lax connections Aj is the gauge transformation A{ = (6±) _1 Af&± + (©±) _1 d{Q±, 
satisfying O-0Z 1 = T^ 1 gTo, where g is an arbitrary constant group element. Assuming that 
©Z 1 = e p (~ l 'e p (~ 2 ' • • • , 0Z 1 = e q (°^ e q<yl ^ e q ^ • • • , where p(—i) and q(i) are linear combinations 
of grade (—i) and respectively, the zero grade component of leads to the solution 

e ff(o) = B- x e~ vC , B = exp . 

From this we have 

<A , | J B- 1 |A> e ^ = <A , |T - 1 5 To|A>, 

where |A) and (A'| are annihilated by the grade (— i) and (i) generators. Taking the highest 
weig ht states |Aj), i = 0, 1 of sZ(2)W we get the following tau-functions 

Tl = e-+- v = (AilTo-^TolAi), r = e~ v = (X \T^gT \X ). 

The so called solitonic specialization corresponds to the situation when g is given by the 
exponential of an eigenvalue F(z), z € C of the operator ad£^ fc ). In this case we have 

[E^,F(z)] =-2z n F(z), 

where F(z) is the vertex operator 



F(z) 



oo 



=— oo 



With the following result 

T 1 gT = exp{p(z)F{z)}, p( 7 ) = exp 2 {-t^z~ 3 - x~ z~ x + x + z + t +3 z 3 } , (39) 
we get the 1-soliton solution depending on the first four times of the hierarchy 

0(i_ 3 , x~, x + , t+s) = - log ( ^ , (40) 

where t\ = l+|p(z), tq = 1— |p(z) and 6 = const. The field (fl0"|) is a simultaneous solution of the 
bosonic limits (with \i = 1) of the equations ([32]) . ([34"[) . ([36]) and the whole set of equations ([37]) 
and (|38l) . The interesting point is that the second and third terms of the RHS of (|37[) cancel 
each other and the field (f> has to obey 

d +3 <9_30 = -2sinh[20], 

which is the case. Then, the solution ()40|) behaves under ^3 like a sinh-Gordon soliton solution 
but with a cubic rapidity. For higher grade times the extension in (I39p is direct. At this point 
we can notice that each chirality of the extended hierarchy is attached separately to the poles 
z = +00 (positive flows) and z = (negative flows) of the Lax pair with spectral parameter z, 
see for instance [7]. This is exactly the pole structure entering the definition of the sinh-Gordon 
Lax pair. 

We end by giving the fermionic currents for the times £±3 and D^ +l / 2 \ After a lengthly but 
straightforward computation, the equation ([26]) becomes 9+3 Jztj 1 ^ + <9_ 3 J, "3 = with the 
current components 



Ij(+ 1/2) = 3-03- Y>- - 4>-d 2 _ cf>, 



u +3 

-^ji + 3 1/2) = - 2 ^+ sinl#] - (l^+d+tf, - (5 + 0) 3 - \d + ^dl4> 



18 



D.M. Schmidtt 



3.2 The N = (2, 2) Landau-Ginzburg Toda model 

J [3] 



Take the psl(2, 2)f?, superalgebra (see Appendix [A] for details). The Lagrangian is given by 



k 

I = -— [<9+0i<9_0i + <9+0 3 <9_0 3 + i^id+ipi + i) 3 d+ip3 + ^2d-ip 2 + ^id-^i ~ V] , 
V = 2 cosh[20i] — 2 cos [203] — 4 (^iV^ + ^3^4) cosh[0i] cos[0 3 ] 
+ 4 (-01^4 - ^3^2) sinh[0i] sin[0 3 ]. 

Taking D(+V2) = jD (+V2) + D (+i/a) with ^C+i/a) = ^+1/9 ^(+1/2) = e2jF ( + i/2) and 

£,(-1/2) = D (-l/2) +jD (-l/2) with ^(-1/2) = _ iF; (-l/2) )jD (-l/2) = - e2F (~m : the supersymmetry 



transformations (fT7|) and (fTSj) with 0>j, ^ — >• > |"0i are 



5+1/201 = +2 (ei-04 - £2^2) > 5 +1/2 3 = -2 (e2"04 + ^1^2) ) 
5+1/2^2 = + (e2<9+0i + ei<9+0 3 ) , 5 +1/2 '04 = ~ (ei<9+0i - £29+03) , 
&f 1/2^1 = -2ei cosh[0i] sin[0 3 ] - 2e 2 sinh[0i] cos[0 3 ], 
5+1/2^3 = +2ei sinh[0i] cos[0 3 ] - 2e 2 cosh[0i] sin[0 3 ], 

(5-1/201 = +2 (ei*0 3 - 62-01) , 5„ 1 / 2 3 = +2 (fl^l + ^3) j 

S-i/2^i = + (e2<9_0i - ei9_0 3 ) , ^-i/ 2 ^3 = - (ei<9-0i + e 2 9_0 3 ) , 

= -2ei cosh[0i] sin[0 3 ] + 2e 2 sinh[0i] cos [03], 
5-1/2^4 = _ 2ei sinh[0i] cos[03] - 2e 2 cosh[0i] sin[0 3 ]. (41) 

We can check (|20p by applying (|41|) twice giving 

[5+i/ 2 , 5 +1/2 ] = + 4 ( e i e 'i + e 2£ 2 ) 9+, [5-1/2, <5li/ 2 ] = +4 (eiei + e 2 e 2 ) d_, 
[5+i/2> 5 -i/2] = 0. 

We have four real supercharges N = (2,2) because dim/C^ 1 ^ = 2. They are extracted 



from (J29J) and (|30p and are given by 
flow 5+i/ 2 : 

= Q\Ft 1/2) + Q 2 + Ft 1/2 \ 



Q+ = j dxl (^2<9+03 - "04^+01 - 2^1 cosh[0i] sin[0 3 ] + 2^ 3 sinh[0i] cos[0 3 ]) , 
Q+= dx l (-0 4 <9 + 3 + ip 2 d+<t>i ~ 2i>i sinh[0i] cos[0 3 ] - 2^ 3 cosh[0i] sin[0 3 ]) , 



flow 5_!/ 2 : 

Q(+l/2) =Q l 1 r i (+V2) + Q2 j p 4 (+V2) j 

Q l _ = J dx l ('0i9_0 3 + ■03^-01 + 2-02 cosh[0i] sin[0 3 ] + 2-0 4 sinh[0i] cos[0 3 ]) , 
GPl = J dx l ('0 3 9_0 3 - ■0i9_0i - 2-02 sinh[0i] cos[0 3 ] + 2^ 4 cosh[0i] sin[0 3 ]) . 
Now introduce the complex fields 

■0- = -01 + ^03, lp+ = ~1p2 + #4) = 03 + »01 
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and the superpotential W ((f)) = 2fj,cos<p in order to write the lagrangian in a more familiar 
form. Then, we have 

V = \W (0)| 2 - [W" (</>) + W*" ((j)*) i>*_i>* + ] , (42) 

where \ W (4>)\ 2 = 2/j, 2 cosh[20i] — 2fi 2 cos[203]. This Lagrangian is invariant under the common 
N = (2, 2) superspace transformations for a complex chiral bosonic superfield. In terms of the 
new complex fields we have for (|4ip with e_ = — (ei + ie%) and e+ = e\ — iZ 2 that 

<^+i/2<A = ~2e-ip+, S-i/2(t> = +2e + tp-, 

Ki/2^+ = +etd + <t>, 5_ 1/2 iP + = -ie + W>* (<f) , 

<5 +1/2 ^_ = -ie.W* (<f ) , 6_ 1/2 ip. = -eld.cf) 



plus their complex conjugates. Define now the following complex combinations of the super- 
charges \Q± = =F<3± + iQ 2 ± 

Q± = 2 J dx l (V±9±0* T i^W* (<j)*)) , Q± = 2jdx 1 (ip*±d±<l> ± iip^W (</>)) . 

Now, with the Dirac brackets {4>,d t <p*} = {<P*,d t (j)} = +1 and {ip+, = {ip-,ipl} = +1/2 
we have with Q±f = {Q±, /} 

= -2ij, + , Q + (j> = 0, 

Q+*P+ = o, Q+i>+ = d+(f), 

Q+i>. = -iW* (</>*) , Q+V- = o, 

Q_0=-2V_, Q_0 = O, 

= iW* (</>*) , = 0, 

Q-i/>- = 0, Q_^_ = d^<p 

plus their complex conjugates. 

Finally, the total variation becomes 5susy = £-Q+ — £+Q- + £*~Q+ — with 

<5susy</> = +2e+^- - 2e_V+, <5 SU sy<A* = +2e* + - 2e*_i/>* +1 

5susyV+ = +e^d + <f> - ie+W* (<f) , Ssvsy^* + = +e.d + 4>* + ie*+ W (</>) , 
JsusyV- = -e+d-<f> - ie-W* (<f ) , ^susyV- = -e+d-<t>* + ie*_W (<j>) , 

which are the usual N = (2, 2) supersymmetry transformations obtained by using a superspace 
approach. As in the case of the N = (1, 1) model, we expect the existence of higher non-local 
fermionic flows for this model. 



Remark 4. The action (|42[) is a Landau-Ginzburg model on a flat non-compact trivial Calabi- 
Yau manifold X, i.e. X = C. As is well known this model is l?-twistable. It would be interesting 
to see the relation of its chiral ring and the chiral ring of a topologically twisted version of 
a superstring on AdS 2 x S 2 where the action (|42p is extracted as a Pohlmeyer reduction [1TJ . 
This would be in principle, a simple way to test up to what point the Pohlmeyer reduction 
can be understood as an equivalence of quantum field theories, at least at the level of ground 
states and by eliminating conformal anomaly issues. It would be also interesting to trace back 
(if possible) the role played by the extended conformal symmetries of W-type in terms of sigma 
model variables. 
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Remark 5. In the relativistic sector of the AKNS hierarchy which is associated to the Ho- 
mogeneous gradation, the Toda potential has a symmetry, i.e. ^ 0. In this case further 
reduction of the model can be performed, eliminating the flat directions. This is done by cou- 
pling to a quadratic gauge field A and is equivalent to the introduction of singular metrics 
defining the non-Abelian Toda models also known as singular Toda models. The minimal cou- 
pling of the fermions with the gauge field gives two terms proportional to . After integration 
of A, the action will have a potential term which is quartic in the fermions, roughly of the form 

which is possibly related to the curvature R of the background metric. This hierarchy is relevant 
from the point of view of Pohlmeyer reductions where reductions of non-linear sigma models 
inside this hierarchy are quite common |17j . 



4 Concluding remarks 

By coupling two identical supersymmetric integrable hierarchies we have shown that the usual 
notion of superspace/supersymmetry is embedded and alternatively described by the symmetry 
algebra spanned by the subset of flows £-i/2 ; t+x/2-, C (i±i/2j *±ij £±3/2; • • • ) ■ We have 
given the explicit form of the supercharges generating these extended (N, N) supersymmetry 
flows and also shown that the higher grade fermionic flows are inevitably non-local on both chiral 
sectors. In particular, when the t_i/2> *+i/2j flows are supplemented by the algebraic 

conditions = 0, the integrable model is restricted to a reduced group manifold spanned 
by the invariant subalgebra of a reductive automorphism r rc d, see Appendix [A] for an explicit 
example. We do not supersymmetrize the fields, i.e. the 'angles' that parametrize the group 
elements G as usually done in the literature. The reduction by r re d should be, in the general 
case, a natural extension to superalgebras of the automorphism used to define the bosonic affine 
(Abelian) Toda models, see [2] for an example applied to the Lie algebra g = sl{m + 1, C). The 
reduction provides a well defined connection between the dressing elements and the physical 
degrees of freedom cf. equation (|lip . as well as the right number of terms in the potential 
appearing in the action functional by truncating it at the second term, thus the natural affine 
superalgebras involved in the supersymmetric affine (Abelian) Toda models are the twisted 
ones. 

What remains to be done in general terms is a formal proof of the statement that the in- 
troduction of r re d is responsible for the locality (Q± = 0) of the lowest supersymmetry flows 
(£±1/2)- This is equivalent to an explicit construction of a reductive automorphism r rcc j in which 

the invariant subalgebra g re( j has no bosonic kernel of grade zero, i.e. /4 0) = 0. Such 

7"red will 

define in principle, all supersymmetric affine Abelian Toda models attached to the mKdV hier- 
archy (in the A series) and it would be interesting to construct it also for other series of affine 
superalgebras. The next step is to introduce a super tau-function formulation which we expect 
to be a natural generalization of the one introduced in |21] for the bosonic affine Toda models, 
where an infinite number of conserved charges were written in terms of the boundary va- 
lues B of a single tau function r, in the form ~ <9± n logT|#. As we are using a fermionic 
version of the Toda models coupled to matter fields constructed in [13], it is possible that one 
has to consider a 'matrix' of tau functions r mn when solving the whole system (I12p . In this case, 
the interpretation of the single r function as a classical limit of the partition function of some 
quantum integrable system will change or will have to be modified in an appropriate way. The 
point is that quantization can be done, in principle, by quantizing T mn , i.e. by using a quantum 
group of dressing transformations. This obviously deserves a separate study. 
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A potential application for our supersymmetric affine Toda models is related to the Pohlmeyer 
reduction of supersymmetric sigma models. In the present construction we needed to impose the 
conditions = in order to have Kg = 0. This means that the corresponding supersymmet- 
ric reduced models belong to the mKdV hierarchy and that no gauge symmetries are involved, 
thus having Toda models of Abelian type (see equation (j4"2j) for an example) . To apply these re- 
sults in the most general situation we need to find a way of coupling two supersymmetric AKNS 
hierarchies in which K,^ ^ 0, i.e. the reduced models will have gauge symmetries, thus being 
Toda models of non-Abelian type. We expect to deduce the action functional for the supersym- 
metric non-Abelian affine Toda models similar to the one constructed in [17} equation (6.49)] 
in the particular case of the reduction of the AdSs x S* 5 superstring coset sigma model. We also 
expect to deduce the supersymmetry transformations by treating them as fermionic symmetry 
flows in the AKNS hierarchy in the same way as it was done for the action f|28|) and (jXTJ) , (|18|) 
in the mKdV hierarchy. It is worth to compare our supersymmetry transformations with the 
ones proposed in [T71 equations (7.21), (7.22)]. We also comment that our Lax pair and action 
functional naturally includes a spectral parameter and that it is also conformal invariant as a 
consequence of the two-loop nature of the affine algebras used. The study of the 'off shell' 
supersymmetric AKNS hierarchy and its symmetries is already under investigation |23| . 

A Used superalgebras solving Q± = 

We consider the reductive automorphism r r0( j for sl(2, 1)W only [3]. It is defined by 

r rcd (E^) = E { ^ VeJ for a € bosonic root, 

r rcd (£W) = i(-l) q \ E ^) E^^ for a € fermionic root, 

where q(E^) is the grade of E^\ i.e. [Q, ] = q(E^)E^\ r]E a is defined by [a ■ H, E a ] = 
ilE^Ea and a are the roots of the underlying finite-dimensional Lie superalgebra. 

Invariance under r re< j define a twisted superalgebra s/(2,l)( 2 ) C sl(2, 1)W. This twisted 
superalgebra solve the conditions Q± = 0. In [15], two sub-superalgebras of sl(2, l)^ 1 ) solving 
the conditions Q±* = where found by another method. The first algebra solution was denoted 

{ 2) 

by s/(2, l)^/ and coincides with the subalgebra invariant under r re d above. The second algebra 

(2) 

solution was denoted by s Z ( 2 , 1 ) [ 2 ] These two subalgebras gives rise to integrable models of 
sinh-Gordon and sine-Gordon type which are related by analytic continuation in the fields. 
Similarly, in the case of psl{2,2) <yl \ four sub-superalgebras solving the conditions Q± = were 

(2) 

found, they were denoted as psl(2, 2)L ; , i = 1,2,3,4 and give rise to integrable field theories 
coupling models of (sinh-sinh), (sine, sine), (sinh,sine) and (sine, sinh) Gordon type respectively, 
all of them related by analytic continuation. The subscript [i] is just a label used to denote 
a particular subalgebra solving the locallity conditions = 0. The reason of using these 
subalgebras is to turn local the lowest supersymmetry flows. The explicit difference, for the 
case sl(2, 1)W (non-local) and sl(2, 1)?? (local) was worked out in detail in [3]. 

In this paper we consider, for the sake of simplicitly and with the aim of not being repetitive, 

(2) (2) 

only the superalgebras sl{2, 1)^ and psl(2,2)^ whose definitions are as follows 

Q = 2d+ 1 -H 1 , Eg* = \2 ■ _ + flgi)) and 

K B = {K? n+1) = A 2 • H^ n+1 / 2 \ Af n+1) = E%> + E^}, 
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M B = {M[ 2n) = H[ n \ M 2 (2n+1) = Eg> - E^}, 



Mi 



F (2n+3/2) = /g. 



+ E (n+l/2)^ + ^(n+1/2) 



+ E (n+1) \ 
' — ai— 02/ 



n (2n+l/2) 



(E 



(2n+l/2) 



(n+1/2) 
012 

n(n+l/2) 



01+02 
(n) 



+ E 



Q-1+CX2) 



fp(") , p (n+l/2)> 
^-02 ^ ^-01-02; 



(2n+3/2) 



>) - (E 



i^aj ^Qi+^y ^ \ -«2 ■ E/ -ai-a2/' 

rp("+l) p(»+V2)\ / p (n+l/2) F (n+1) 
V^«2 - ^01+02 ) ~ \ h -a 2 



E (n+1) \ 
— ai— 02/ 



S Z(2,2){ 2) : 



Q = d+-(H 1 + H 3 ), 



n(±l) 



and 



Mb = 



M 



(2n) 



r(2n) A/f (2n+l) 



p (2n) _ p (2n+2) 
"Qi - C/ — ai i 



M (2n + l) =iJ (2n + l) ) M (2n) 



(2n+l) 



F 



(2n+l/2) 



Ea3 

(2re+l/2) 



/ p (2n+l/2) p (2n+l/2) n / p (2,«.- v ^ 
^a 2 "I" -C/ Ql + Q2 + a3 ; "I" ^-C/q-j^^v^ 



,(2ra-l) 
03 

,(271-1/2) 



/ p (2n+l/2) „ 
' > — Q2 — ai — a.2— ocz 



? {2n+\/2) 



1 -^02+03 ) ~'~ 
(2n+3/2) „(2n-l/2)> 



,(2n+3/2)> 
02+03 



\ /p(2n+3/2) p (2n-l/2)N 
/ ' \ —ai—aa ' — 012— 03 /' 



p (2n+3/2) _ / p (2n.+3/2) p (2n+3/2) s / p (2n+l/2) p (2n+5/2)x 
-To — 1-C/Q2 I" -C/ ai+a , 2+0 , 3 ; "I" y& ai +a 2 "T ^02+03 / 



/ F (2n+3/2) F (2n+3/2) n _ / p (2n+5/2) „(2n+l/2h 
V — Q2 ^-01—02-03/ V — Q!l— «2 — "2— «3 /' 



(2n+3/2) 



/ p (2n+3/2) p (2n+3/2) v _ / „(2: 



(2n+l/2) . ^(271+5/2) 



\-OL2 



'02+^3 



») + 



+ (E 



-a 2 



F 



(2n+l/2) 



+ E[ 
/„(2n+l/2) 



,(2n+3/2) _,_ F (2n+3/2) ^ _ /j ? (2n+5/2) „(2n+l/2)\ 

I — ai — a 2 —0.2—0.3, ) ' 

P^jh-i/zj \ /„(2n-l/2) p(2n+3/2)\ 

«l+02+Q!3/ V Q1+CK2 "T" - C/ 2 +Q3 / ~ 1 ~ 



TWf = < 



G 



+ ( E 

(2n+l/2) 



,(2n+l/2) 



-a 2 



+ E 



(2ra+l/2) 

,(2n+l/2) 
-a% — 02—03 



s _ /p,(2n+3/2) F (2n-1/2)N 
/ V — 01—02 ~ l ~ — Q2— «3 / 



,(2n-l/2) _ p (2n+3/2)N / p (2n+l/2) 



<2n+3/2) 
-ai— a?2 



(2n-l/2) 



p (2n+l/2) \ 



-02— "3 



(2n+l/2) 
«2 



^(2n+l/2) ^ 



-ai— «2— «3/ 



^(2n+3/2) _ ^(2n+l/2) 



J a\+a2 
/p(2n+5/2) _ 
\ — ai— a2 



,(2n+5/2) 
02+^3 



)+(E 



,(2n+3/2) 
'02 



(2n+3/2) 



(2n+l/2) 
-Ct2~ «3 



(2n+3/2) _ £,(2n+3/2) 



-a 2 



a.\—ot2—ot'ij 



G 



(2n+3/2) _ /£,(2ra+l/2) F (2n+5/2) 



+ (E 



G 



,(2n+l/2) 



ai+a 2 
,(2n+5/2) _ 

Cti — 02 

,(2n-l/2) 
ai+«2 



Z A / P (2n+3/2) p (2n+3/2) \ 



E 



'02+03 

(2n+l/2)\ _ / p (2n+3/2) 
-02-03 / V -02 

,(2n+l/2)^ , ( F (2»i+l/2) 

'Q2 



(2n+3/2) 
01—02—03 



r (tn-i/zj j ? (2n+l/2)\ / p 



■ E 



). 

,(2n+l/2) 
01+02+03 



) + 



(2n+3/2) 
-ai— 02 



(2n- l/2)i 
-a 2 — a 3 , 



(2n+l/2) ^(2n+l/2) 



01— Q2— 03 



3) 



B Relativistic equations for t-1-3 

The decomposition of the zero curvature equation [L + 3,L_3] = on the different graded sub- 
spaces reads 



5±3jD (+5/2) = ± ^( T3))i? +l^(±l/2) j g±l] i 

d± 3 D^) = ±[dW 2 \b* 1 d( ±1 Mb ±1 ) ± [E^^Ti^Cii)^^ 

5±3Z) (+3/2) = ± ^(+2) )BT ljr,(±l/2) 5 ±l] ± [^(TS/2) ^+^(±1)^+1] 
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± [ £ ;(T3) jB Tl £ ,(±3/2) B ±lj j 

5 ±3J D (Tl) = ±[dW 2 \B* 1 D( ±1 MB ±1 ] ± [£(^^^£,(±1)^1] 
± [£,(T5/2) jS Tl£,(±3/2) B ±l] ± [£(^£^£,(±2)^1^ 

± [D^ £^£,(±3/2)^1] ± [D^s/^^i^i^ii] 
3„ 3 (d+sBB- 1 ) = [£)(+ 1 /2) )jB£ )(-V2) jB -i] + [£,(+i) j5£) (-i) B -i] 

where we have replaced the solution = —d+^BB -1 for the degree —3 equation 

(43) 



D(-V2) =V ,l G f(-6/2) j 

d(-2) = 4m 1 ( ~ 2) , 

£,(-3/2) = ^^(-3/2) + ^3 G (-3/ 2)) 

= ^^-^ + 3 K^ l) + ^M 2 ( - 1} , 

D(-V2) = ^4 jP (-l/2) +V; 5 G (-l/2) ! 



5±30^ = 2^V± cosh[</>] ± 244 sinh[0] - 24 sinh[20] ± 24 cosh[20], 
3±s4 = ~44 smh[(j)} t 44 cosh[0] + 44 sinh[20] t 44: cosh [20], 
d± 3 t/4 = =F (44 - 2^|) cosl#] - (44 - 2</4) sinh[0] ± 44 cosh[20] 

-44sinh[20]±44, 
d±s4 = 2 (~44 + 4^± - 44) cosh M ± 2 ("44 + 44 " 44) sinh[0], 

£>± 3 4 = 2 (-44 - 44 + 44) coshM t 2 (44 + 4^± - 44) 

- 244 sinh[20] ± 244 cosh[20], 

a ±3 4 = ±2 (44 + 44 - 44) sinl#] + 2 (4 4 + 44 - 44) cosh[0] 

± 244 cosh [20] - 244 sinh[20] ± 24, 
^4 = =F (44 " 44 - 44 + 44) cosh[0] 

+ (44 - 44 + 44 - 44) S mi#] 

t (44 + 44) cosh t 2 ^ + (44 + 44) ^[^] ± 44 ± 44, 
5 ±3 4 = t (44 + 44 + 44 + 44 - 24) cosi#] 

- (44 + 44 + 44 + 44) S mi#] - (44 + 44) smh[20] 
± (44 + 44) coshp^] ± 44, 



= [D (0) + d+sBB' 1 , BE^B- 1 } . 

(2) 

For the s/(2, l)^-/ algebra set 

£)(+ 5 /2) = ^ G (+5/2) ) 

l>(+ 2 ) = 4m{+ 2 ), 

D(+3/2) =V ,2 F ( + 3/2) +V; 3 G ( + 3/ 2 ) ) 

= 4^i +1) + 4^ +1) + 4m 2 (+1) , 

D(+V2) = ^4 i?2 (+l/2) +i/; 5 G ,( + l/2) ) 

z>(°) = 4m} ) 

then, the above equations of motion are 
<9 ±3 4 = 24 sinh[0] ± 24 cosh[0], 
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+ 2 (i/>%il)i + i/jiipl - V^^l - cosh[0] 

- 2 (1 - + 0^1) sinh[2$ - 2 ($.01 _ ^4 ) cosh [ 2 $. (44) 

The generalized equations of motion are found by replacing above the following solutions 

4 = -\d%<p + i (a+0) 3 + zd + ^ + d + ^ + . 
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